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ABSTRACT 

We use the structure of the Einstein ring image of the quasar host galaxy in the four-image quasar lens 
PG1 1 15+080 to determine the angular structure of the gravitational potential of the lens galaxy. We find that it 
is well described as an ellipsoid and that the best fit non-ellipsoidal models are consistent with the ellipsoidal 
model. We find upper limits on the standard parameters for the deviation from an ellipse of |af /a$\ < 0.035 
and |«4 /oo| < 0.064. We also find that the position of the center of mass is consistent with the center of light, 
with an upper limit of | Ax\ < 0."005 on the offset between them. Neither the ellipsoidal nor the non-ellipsoidal 
models can reproduce the observed image flux ratios while simultaneously maintaining a reasonable fit to the 
Einstein ring, so the anomalous flux ratio of the Aj and A2 quasar images must be due to substructure in the 
gravitational potential such as compact satellite galaxies or stellar microlenses rather than odd angular structure 
in the lens galaxy. 

Subject headings: cosmology: observations — dark matter — gravitational lensing — quasars: individual (PG 
1115+080) 



1. INTRODUCTION 

The scenario of hierarchical structure formation dominated 
by cold dark matter has been extremely successful in explain- 
ing various obse r vations of large scale structure (see, e.g ., 
lYork et alJ 120001: iNorberg etalTl2001t ISoergel et all 120031) . 
In this scenario, small structures f orm first and merge to 
build progressively larger structure llRees & Ostrikerjfl977t 
IWhite & Reesl [l978). One prediction of these theories, both 
analytically and in simulations, is the existence of far more 
low-mass halos than are observed either in the field or as 
satellites of the Milky Way. Detailed studies of small satel- 
lite halos show that many survive without being tidally dis- 
rupted, albeit with some mass loss, and remain as gravita- 
tionall y bound substructures in the ha lo of the larger galaxy 
JMoore et al]ll999tlKlvpin et alJI 1999ft . Most estimates lead 
to an abundance of satellites that are inconsistent with the ob- 
served abundances, suggesting that the substructure must con- 
sist of dark satellites in which star formation was suppressed 
by heating and ejecting the baryons ( Bul lock et alJ2000T) . Un- 
fortunately, dark satellite halos are difficult to detect. 

Fortunately, we have one probe of halo structure that can 
detect substructure purely thro ugh its gravitational effects - 
gravitational lensing (see, e.g., iBl andford & Naravanl ll992t 
ISchneider et all 19921 iKochanek et all2004l) . Substructure in 
the halo of a gravitational lens galaxy, whether dark satellites 
or simply stars, can modify the observed fluxes of the individ- 
ual images relative to those in a smooth gravitational poten- 
tial. These are most striking for the systems with anomalous 
flux ratios between images which should have similar flux ra- 
tios given an y smooth potential center e d on t he prima ry lens 
galaxy (e.g . ISchneider & Maol IZhao & Pronkl 120011: 

Chiba 2002; Dalai & Kochane kl2002HMetcalf & Zhaol2002l) . 
For the case of microlensing by the stars in the lens galaxy, 
the existence of substructure is easily verified by the tem- 

Based on Observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated 
by AURA, Inc., under NASA contract NAS5-26555. 



poral variatio ns of the flux ratio s uch as those observed in 
Q2237+0305 JWozniak et alJ EOQO). Tests for satellites have 
been focused on lensed radio sources that are too large to be 
affected by microlensing and little affected by th e interstellar 
medium of the lens galaxy. Kochanek & Dalai (2004) have 
shown statistically at least, that the flux ratios of lenses show 
the characteristic pattern of demagnified saddle point images 
expected for low optical depths of substructure embedded 
in a higher optical depth of smoothly d istributed dark mat- 
ter ( ISchechter & Wambsganssll2002l lKeetonll2003l) . If these 
anomalies are interpreted as being due to s ubstructure, the es- 
timate d mass fraction in the substructure ( Dalai & Kochanek 
2002) is consistent w ith theoretical expectations ( Moo re et alJ 
1999; Bullock et al. 2000). There is some debate over the 
contribution from small halos along the line of sight relative 
to those associated with the lens galaxy, with most studies 
favoring a dom inant contribution from the lens galaxy (e.g., 
lMetcali2002llChen et alJl2003l) . However, from the point of 
view of understanding the halo mass function this is a moot 
point because these low mass halos are not observed as either 
satellites or in the field with an abundance approaching that 
predicted by the CDM scenario. 

There have, however, been several explorations of whether 
the flux ratio anomalies can be explained by adding complex 
angular structures to the gravitational potential of the pr imary 
lens galaxy rather than by adding substructure llEvans & Witll 
120031 Kochanek & Dalai [2004 iMoller et alJl2002l) . In pa7 
ticular, Evans & Witt (2003) showed that for one particularly 
simple model with arbitrary angular structure the fit to the 
data can be reduced to a problem in linear algebra in which 
the lens constraints, including the flux ratios, can be modeled 
to arbitrary accuracy simply by including enough terms. The 
model parameters of this study were problematic as the quasar 
images needed only to be fitted to 0"05, which is more than 
ten times the actual uncertainties. However, the mathemat- 
ical observation that standard models fitting only the image 
positions fail to fit the image fluxes only because the gravita- 
tional potentials are assumed to be roughly ellipsoidal is cor- 
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FIG. 1 . — Einstein rings extracted from the data subsamples described in 
O Heavy solid and dotted lines represent the radius of the Einstein ring of 
the final combined data and its uncertainties, respectively. The four image 
positions are shown as filled circles. 



rect. However, iKochanek & Dalall d200l considered models 
with some additional angular structure and compared models 
constrained by either only the positions of the quasar images 
or the positions of the images combined with any additional 
constraints available for the system. In each case, the large 
deviations from an ellipsoidal model that could improve the 
modeling of the flux ratios were ruled out by the additional 
constraints. 

In this paper we explore whether the lens galaxy in the 
four-image quasar l ens PG1 1 15+080 (IWevmann et alJU98?)t 
Courb iyetatlli997t iKeeton & Kochaneklll997t llmpev et all 
1 19981 iTreu & Koopmans 2002) has significant deviations 
from an ellipsoidal density distribution. PG1 115+080 is an 
interesting case because the merging A1/A2 images have a 
flux ratio of ~ 1 .6 rather than the ratio of ~ 1 .0 predicted by 
smooth models. We will use the same models with arbitrary 
angula r structure considered by Evans & Witt (2003, and ear- 
lier by IWitt et alJl2000t lEyans & Wit3l200U IKochanek et al.1 
120011 IZhao & Pronkl l2001) but examine the constraints im- 
posed by the Einstei n ring im age of the quasar host galaxy 
in detail. IKochanek et al.l (1200 ll) demonstrated that the shape 
of the Einstein ring completely encodes the angular structure 
of the potential, so we can quantitatively measure the devia- 
tions of the higher order multipoles of the lens gravity from 
the best fit ellipsoid. In 5j2]we present the Hubble Space Tele- 
scope (HST) observations on which our models are based, and 
then briefly review our models and the fitting procedures in 
|3] Our main results are presented in 50 an d we summarize 
them in |j5] 

2. HST OBSERVATIONS 

We use the two HST //-band (NICMOS/NIC2/F160W) im- 
ages of PG1 1 15+080. The first observation was a single-orbit 
(2 560s) observation tak en on 1997 November 17 and reported 
bv llmpev et al.l d!998h . The second observation was a four- 



orbit observation (11700s) taken on 2003 June 15. For both 
observations, we used four dithered sub-images per orbit to 
remove hot pixels and to imp rove the flat-fie ldin g. The data 
were r educed as described in lMcLeodl i 19971) and lLehar et alJ 
(2000). The spacecraft orientation was 69° for the first ob- 
servation and -112° in the second observation, so there was 
little change in the orientation of the quasar diffraction spikes 
relative to the lens. 

We e xtracted the ring curves as described in Kochane k et alJ 
(HooB. We first fit a photometric model to the image includ- 
ing the four lensed quasar images, the lens galaxy and the host 
galaxy. This fit provides estimates of the quasar image fluxes, 
and these are used to subtract the lensed quasar images from 
the images. We then select a point close to the center of the 
lens, and as a function of position angle 9 we find the radius 
r(9) of the peak surface brightness of the Einstein ring and 
its uncertainties. The total number of points extracted from 
an Einstein ring is A^i ng = 71. Since the ring is well-detected 
even in a single image, we extracted rings from several sub- 
sets of the data to explore both statistical and systematic er- 
rors in the data. We will refer to extractions from the 1997 
data as obsl and extractions from the 2003 data as obs2. For 
obsl we extracted the ring using three different choices of the 
point about which to extract: obs 1 - A was centered on the lens, 
obsl-B was offset by 0."1 in the x direction of the image and 
obsl-C was offset by 0"1 in the y direction. For the obs2 data 
we extracted the ring for each of the four individual orbits 
(obs2-A, obs2-B, obs2-C and obs2-D). For our final analy- 
sis we used a ring extracted from the combined image of the 
four individual orbits in obs2. Figure^shows r{6) for these 
various extractions as compared to the final combined model 
and its uncertainties. For the image positions an d fluxes, we 
use the //S/7NICMOS observations described by Impe v et all 
(1998). 

3. LENS MODEL OF PG1 115+080 

In this section, we first describe our models for the lens 
galaxy of PG1 1 15+080 and the nearby galaxy group as an ex- 
ternal perturber. Then we describe our model for the Einstein 
ring and our approach to fitting the data. 

3.1. Lens Potential 

We model the lens galaxy as a scale-free potential with ar- 
bitrary angular structure, 

f(0) 



<Kr,0) = rF(0), K (r,9) = 



2r 
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where F(9 ) is an arbitrary fun c tion of angle 9, and f(9) 
I + F" dZhao & Pronkl l200ll lEvans & Wittl I2OO1I 12001 
IKochanek et al.ll2001t IWucknitzl 120021) . These models have 
a flat rotation curve, making them consist ent with most esti- 
mates of the radial structure of len ses (e.g. JTreu & Koopmans! 
2002; Rusin & Kochanekl 120041) . We expand the angular 
structure in a multipole series, 
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F(9) = flo + [ fl m cos m @ + b m sin m9] , 
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neglecting the dipole (m = 1) terms which will be degenerate 
with a shift in the source position. From Poisson's equation, 
the angular structure of the convergence k is related to that of 
the potential by 
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/(6') = flo + ^^ [a m (l - m 2 ) cos m9 + b m (l — m 2 )sinm#] 
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1 + (A m cos m6 + B m sin m6) 
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where A m = (l—m 2 )a m /ao and fi„, = (l-m 2 )b m /ao for m > 2. 
We are interested in the deviations of the potential from an 
ellipsoidal surface mass density, 

oo 

/(0)oc [l-(l-^ 2 )cos 2 6»]~ 1/2 oc l+J2 a 2m cos2me > ( 4 ) 

where q is the axis ratio of the ellipse with coefficients a\ m to 
distinguish them from the coefficients fl2,„ of the general po- 
tential. We will look at deviations of the gravitational poten- 
tial from an ellipsoid defined by the aq, «2 and Z?2 coefficients. 
In this model, the quadrupole of the ellipsoid is defined by 
3 

a\ = (fl2Cos26'z, + fe2sin2^), (5) 

a 

for axis ratio q and major axis orientation 9 L = 
0.5tan -1 (£>2/fl2)- The remainder of the coefficients a\ m 
for m > 2 of the ellipsoidal model are defined by the standard 
multipole expansion of equation (0}. We describe the rest of 
the model by the deviations Aa m and Ab m (m > 3) from the 
ellipsoid. Thus, our overall model for the primary lens is 

F(9) = ao + fl2cos 2d + bi sin26' 

oc 

m=2 
oo 

+ *^^[Aa m cos,m6 + Ab m sinm8] , (6) 

m=3 

where «o, «2, k>2, Aa m and Ab m are model parameters. 
The deviations in the mass density from the ellipsoid are 
AA m = (l-m 2 )Aa m /flo = flm/ao and AB m = (l-m 2 )Ab m /ao = 
b\ja$ for m > 3 where af n and bf n are the standard coeffi- 
cients that describe the deviatio ns of the isophotes of ellipti- 
cal galaxies from el lipsoids (e.g., Bender & M ollenhof fll987t 
Bende7ItIi]E989). Since the lens galaxy is fairly round 
qi ~ 1, we expand the ellipsoid to order m = 5 and usually 
consider deviations for m = 3 to m = 5 beyond which the coef- 
ficients a\ m are negligible and the perturbations are well below 
the noise. 

3.2. The External Perturber 

When modeling gravitational lens systems, it is important 
to include an independent external shear that can be mis- 
aligned with the light distribution (Keeton et alJll997l) . In 
the case of PG1 1 15+080, the dominant tidal perturbations are 
associ ated with the g roup to which the lens galaxy belongs 
(e.g.. lImpevet"aiTT9 98). The group is composed of two lumi- 
nous spirals and several smaller galaxies with a luminosity- 
weighted centroid ~ 20" from the lens. Since the higher or- 
der perturbations of the group beyond an external shear are 
quantitatively importan t, we need to generalize the models of 
iKochanek et al.l fcOOll) to include these higher order terms. 

We model the group of galaxies as a singular isothermal 
sphere (SIS), (f) ext = b g \r-r g \, where b g is the Einstein radius 
of the group and r g is the position vector of the group relative 
to the lens center. If the group is far away from the lens in 
the sense that b/r g <C 1, then we need to include only the 
quadratic terms in an expansion of the potential, 
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The first term provides t he mass sheet degeneracy and af- 
fects only the time delays ( Fa lco et alJl985UGorenstein et al.l 

■ bj2r„~0.1 



1988). The second term is the tidal shear, 
that i s the dominant external perturbation ( Sc hechter et alJ 
(19971) . However, in PG1 1 15+080 the higher order terms are 
significant, so we must use a more complex model. We ap- 
proximate the potential using only linear and quadratic terms 
in radial dependence, 

oc oo 

&ext(r, 0) ~r 2_\. a m cos tn(6- 9 g ) + r 2 b s m cos m(6 - 8 g ) 



m=\ 



m=l 



rF ext (O) + r 2 G ext (0), 



(8) 



because it has the advantage that the structure of the Einstein 
ring can be analytically calculated even when the potential has 
arbitrary angular structure. The quadratic term of the potential 
Gext(^) is added to yield a better approximation to the radial 
deflections of a true SIS model. We determined the coeffi- 
cients af n and bf n for the external perturber by minimizing the 
difference in the deflection of our approximation from that of 
a true SIS near the Einstein ring of PG1 1 15+080. We need to 
expand <fi ext only for m < 3, and the rms deflection difference 
between the resulting model and a true SIS model at r g ~ 10" 
from the lens is only 0."0005 over a 1" annulus encompassing 
the ring. This is significantly less than our smallest astromet- 
ric errors (0."005). 

3.3. Einstein Ring Models 

We model the Eins tein ring using an e xtension of the the- 
ory developed by Kochaneket al. (2001) for our more gen- 
eral potential, 0(r,9) = r[F(9)+F ext (9)] + r 2 G ext (6). When the 
source is extended, the tangentially stretched images of the 
source form an Einstein ring. If we assume that the source 
has a monotonically decreasing surface brightness, the Ein- 
stein ring radius relative to the lens is simply 



h-S t + U{) S t 



(9) 



t-St 

where h = (F + F ext )e r + (F' + F^ xt )ee, the source plane tangent 
vector is t =M~ l ■e r = (l-2G ext )e r -G' exl ee,M~ l is the inverse 
magnification tensor, the source center is «o, and 5 is the two- 
dimen sional shape tensor of the source (see Kochan ek et alJ 
2001] for the details). The only relevant parameters of the 
source are its axis ratio q s and the position angle of its ma- 
jor axis 9 S . 

3.4. Model Fitting 

We use a simple \ 2 statistic for the goodness of fit of a 
model described by parameters p. The astrometric constraints 
consist of the observed Einstein ring radius r(6*,) and the four 
source positions u,. We divide the estimate into the goodness 
of fit to the astrometric constraints \l s t an ^ me f ux constraints 

Xlu 



1 The fit statistic for the astrometric constraints is 



xL(p) = Yl ' 



[u -Ui(p)]Mi\ 



^ \r(ed-r(6r,p)\ 2 
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0~7 



(10) 



/.ring 



where M,- is the magnification of image i, and the cr,'s are the 
uncertainties in either the image positions or the Einstein ring 
radii. The fit statistic for images is approximately correct be- 
cause of the magnification weightings - the need to analyze 
models quickly in a large parameter space precluded using 
the exact image plane fit statistic. We can also constrain the 
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FIG. 2. — A check for systematic errors by comparing the results for the 
rings extracted from the data subsamples listed in the top left panel and dis- 
cussed in For illustration, we show the shear and PA of the group (top 
left), the lens galaxy position (top right), the quadrupole moment of the den- 
sity distribution (lower left) and the 58 moments (lower right). The formal 
uncertainties in parameters are only shown for the combined data. For statis- 
tical uncertainties, the scatter in the subsamples should be twice the errorbar 
for the combined data. 



from the four individual orbits of obs2 (obs2-ABCD). In each 
case the lens was re-modeled as part of the extraction pro- 
cedure. Unfortunately, the rotation of the PSF between the 
two observations is negligible (see |j2j, so comparisons of the 
data sets provide little constraint on systematic errors arising 
from the PSF. Figure |2] compares the results for these seven 
data sets for a range of variables using the full non-ellipsoidal 
models. The main point to note is that the parameter estimates 
are mutually consistent. Also note that the distributions of the 
parameter estimates are consistent with the formal statistical 
uncertainties, except for the 5th order deviations where they 
are smaller. For the remainder of the paper we consider only 
the combined data. 

4. RESULTS 

In this section we address three issues. First, is the lens 
consistent with an ellipsoid when we fit only the astrometric 
constraints? Second, is the center of mass consistent with the 
center of light? Third, are any of these conclusions changed 
if we try to fit the flux ratios as well? 

4.1. Is the Lens an Ellipsoid? 

We start by fitting only the astrometric constraints, first 
with a purely ellipsoidal model, and then adding deviations. 
The results are presented in Tabled An ellipsoidal model 
with a very round lens (qi ~ 1), a group located at a po- 
sition angle of 9 g ~ -113° and at a distance of r g ~ 11" 
provides a reasonable fit to our data and is consistent with 
the earlier models using a singular isothermal ellipsoid for 



model using the image fluxes, although fluxes can be contam- 
inated by systematic errors, dust extinction, microlensing by 
stars, or substructure. For parity-signed fluxes /,-, we add 

xL(WJ) = E^^, (id 

where the source flux f s is an additional free parameter of the 
fit. 

We start by fitting the astrometric constraints with a fiducial 
ellipsoidal model that has 12 free parameters; the source po- 
sition is «o = (u x , w y ), the axis ratio q s and major-axis position 
angle 6 S of the source, the shear amplitude and the position 
of the external perturber {7^, r g , 9 g }, the ellipsoidal lens po- 
tential {oo, 02, i>{\, and the lens center is (x^, y£). Once we 
have estimated the goodness of fit for this model, we can add 
the parameters for the deviations from ellipsoidal symmetry 
or xi U x' an( l examine the change in the goodness of fit and the 
parameters. 

We minimize the x 2 of a given model and estimate the 
parameter uncertainties using Levenberg-Marquardt method 
(e.g.. lPress et alJl992l) . We ensure that we have found a gen- 
uine x 2 minimum by repeating the minimization with a range 
of initial parameters and by checking the existence of the min- 
imum with the downhill simplex method. We tested the code 
on a range of synthetic lenses with Einstein rings including 
both ellipsoidal and non-ellipsoidal models. 

As a last step before presenting our results, we fit rings ex- 
tracted from subsamples of the data as a check for systematic 
errors, and the level of our random errors. We consider the 
standard ring from obsl extracted in three different positions 
relative to the lens galaxy (obsl-ABC) and the rings extracted 
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FIG. 3. — The best fit ellipsoidal (heavy lines) and m = 5 non-ellipsoidal 
(light lines) models for the ring based on the astrometric constraints. The ob- 
served positions of the Einstein ring and the lensed images are represented as 
open circles and triangles, respectively. The line sizes of the open circles are 
the positional uncertainties of the Einstein ring radius while the uncertainties 
in the image positions (0"003) are omitted. The filled circles represent the 
Einstein ring positions of the best fit ellipsoidal model, while the filled ellipse 
at the center shows the predicted position, ellipticity and position angle of the 
source. The heavy solid and dashed lines show the critical line (solid) and the 
isodensity curve (dashed) of the ellipsoidal model, while the corresponding 
curves of the best-fit non-ellipsoidal model are shown with the light lines. 
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TABLE 1 
Best-Fit Models 





Ellipsoidal Model 


Non-Ellipsoidal Model 


Parameter 


Astrometry Only 


Astrometry+Flux 


Astrometry Only 


Astrometry+Flux 


Is 


0.69(02) 


0.69(02) 


0.70(02) 


0.68(03) 


$,(degrees) 


-32.0(2.7) 


-30.9(2.6) 


-33.4(4.6) 


-30.9(3.6) 


7s 


U. 1 1D(U1 j 


U. 1 1Z(U1 j 


U. 1 Izllo ) 


U. 1 1 LyvtZ. ) 


r g (arcsec) 


10.8(1.4) 


12.4(1.8) 


16 9H2 51 


12.5(7.1) 


S ? (degrees) 


-113(00) 


-116(00) 


-113(01) 


-116(01) 


.v/Xarcsec) 


0.001(02) 


0.005(02) 


-0.003(03) 


0.008(03) 


y/,(arcsec) 


0.001(01) 


0.001(01) 


-0.002(03) 


0.004(03) 


c lL 


0.97(02) 


0.98(01) 


0.96(03) 


0.99(02) 


A1/A9 


0.97(03) 


0.94(02) 


1.00(13) 


0.93(06) 


0/,(degrees) 


-66.1(7.8) 


68.4(4.9) 


-52.4(9.2) 


61 8(26 11 






I 1 5 TOO} 


1.15(01) 


1 1 5f001 


10 3 a 2 


4.1(2.0) 


-3.2(0.9) 


2.2(2.7) 


-1.0(1.7) 


10 3 fo 2 


4.5(1.1) 


3.0(0.6) 


8.3(3.4) 


1.5(1.5) 


10 3 Aa 3 


= 


= o 


-1 4(4 6) 


0.4(1.7) 


10 3 A/> 3 


= 


= 


0.8(7.2) 


0.5(2.7) 


10 3 Aa 4 


= 


= 


-0.1(4.6) 


-0.7(1.3) 


10 3 Ai) 4 


= 


= 


0.3(5.7) 


-0.3(1.2) 


10 3 Ao 5 


= 


= 


-0.2(1.9) 


0.4(1.6) 


10 3 Afe 5 


= 


= 


0.3(1.1) 


0.1(0.7) 




35.6 


38.1 


29.3 


41.6 




0.2 


5.4 


0.1 


5.4 




0.1 


2.8 


1.6 


9.2 


Xflux 




45.4 




28.7 


Xtol 


36.0 


91.7 


31.1 


84.8 



NOTE. — The best fit ellipsoidal and non-ellipsoidal (m = 5) models of PG1 1 15+080, fitting 
either only the astrometric constraints or the astrometric constraints and the observed flux ratios. 
All the angles are standard position angles while the lens position and the coefficients of the lens 
potential are calculated in Cartesian coordinates in which x-direction is toward west. The errors 
are shown in parentheses. 



the lens tg chechter et al] 119971 Keeton & K ochanekl [l997 ; 
I Courbin et alJll997| llmpev et aljll998t IZhao & PronklhoOl ; 
kochanek & Dalall 120041) . The total xL - 36 for 69 de- 
grees of freedom means that we have somewhat overfit the 
data and/or overestimated the uncertainties by 38%. Figure[3] 
shows the best fit ellipsoidal model using only the astrometric 
constraints. 

Next, we fit the data using a sequence of non-ellipsoidal 
models to investigate how much the lens can deviate from el- 
lipsoidal. We did this by sequentially including higher order 
poles starting with m = 3 (Aaj, A/73), then m = 4 and finally 
m = 5. For comparison with the ellipsoidal model, we show 
the critical line and the isodensity curve of the best fit m = 5 
non-ellipsoidal model in Figure [3] and we present their final 
parameters in Tabled It is particularly interesting to note that 
all the higher order deviations are consistent with zero. 

We summarize the resulting \ 2 of each model in Table [2] 
We use the F-test to estimate whether the new parameters sig- 
nificantly improve the fit, and we find that none of the new 
variables significantly improves the models. In many cases, 
adding the new variables even raises the \ 2 P er degree of free- 
dom. All the best fit non-ellipsoidal models are consistent 
with the best fit ellipsoidal model. In short, the constraints 
on the angular structure from the Einstein ring force the lens 
of PG1 115+080 to be consistent with an ellipsoidal density. 
The isophotal deviations of the best fit non-ellipsoidal mod- 
els are ~ 10~ 3 for both af n and (m < 5) somewhat smaller 
than the values of af /a Q = 0.009, a% /a = -0.004 for H-bmd, 
and a^/flo = -0.015, a%/a = -0.004 for /-band found by 



Evans & Witt (2003). 

In our first models, we constrained the center of mass to 
agree with the obse rved lens position and its formal uncer- 
tainties of 0"003. Evans & Witt (2003) were concerned that 
the center of light and mass may differ, and allowed the lens 
center enormous freedom to move relative to the observed po- 
sition (0"05). For our next experiment, we drop the constraint 
on the lens position and re-optimize the models. We find 
a lens position of (-0."001 ± 0"002, -0"001 ± 0"002) com- 
pared to the measured position of (0."000 ± 0"003, 0."000 ± 
0"003) or a net difference of 03a for the RA and 0.4cr for 
the Dec. The center of mass of the lens is essentially iden- 
tical to the center of light. The observed lens galaxy of 
PG1 115+080 h as an axis rat i o q ~ 0.9 a nd a position angle 
9 ~ -70° (Imp ev et al]ll998l iFalco et alJl2000l) . The best-fit 
axis ratio qi of the mass is rounder than the light while the po- 
sition angle 61 is aligned with the light for both the elliptical 
and non-elliptical models. 

4.2. The Flux Ratio Anomalies 

Because we believe the anomalous flux ratios are due to 
substructure (satellites or stars) rather than a problem in the 
lens model, we have so far neglected the flux ratios as a 
model constraint. We now investigate the changes in the lens 
potential when the flux constraints are added. The best fit 
parameters of an ellipsoidal model with flux constraints are 
also shown in Table ^ There is little change from the pa- 
rameters of the ellipsoidal model based only on the astromet- 
ric constraints. As expected, most of the contribution to the 
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TABLE 2 

F-Tests of The Additional Degrees of Freedom 



Model (m) 


Astrometry only (A) 


Astrometry+Flux (F) 


N par 


dof X 2 d0f P(%) 


Npar 


dof 


x z dof P(%) 





12 


69 0.522 100 


13 


72 


1.273 100 


3 


14 


67 0.474 69.1 


15 


70 


1.310 90.5 


4 


16 


65 0.488 78.7 


17 


68 


1.301 92.7 


5 


18 


63 0.493 82.0 


19 


66 


1.285 96.6 



NOTE. — Using astrometric constraints only, models with higher order de- 
viations (Am) for order cosmd are compared to the fiducial ellipsoidal model 
(AO). When flux constraints are added, models (Fm) for order cos m(9 are com- 
pared to F0. The dof column gives the number of degrees of freedom and the 
Xjof column gives the x 2 P er degree of freedom. The probability P gives the 
F-test probability that a given non-ellipsoidal model is consistent with ellip- 
soidal model. A non-ellipsoidal model with lower P implies a more significant 
improvement relative to the ellipsoidal model - none of the improvements is 
significant. 



fit statistic Xnux comes from the merging A1/A2 image pair 
showing the well-known flux anomaly. To reproduce the ob- 
served flux ratios, the critical line should be either distorted 
or moved closer to the brightest image so that the flux ratio 
becomes higher than predicted by ellipsoidal models. In the 
lEvans & Wittl (|2003) models, this was possible because the 
quasar and lens positions could be shifted by far more than 
their actual uncertainties. However, with the correct QSO un- 
certainties and the Einstein ring constraints, especially on the 
angular structure of the lens potential, the best fit model with 
the flux constraints largely coincides with that using the astro- 
metric constraints only. 

We then added the higher order deviations from an ellipsoid 
and fit the data with non-ellipsoidal models. Since the mag- 
nification, M oc F" oc m 2 a m , is more sensitive to the higher 
order multipoles, the best-fit non-ellipsoidal models produce 
a better fit to the flux ratios, but it is still impossible to signif- 
icantly improve flux ratios while simultaneously maintaining 
a good fit to the Einstein ring. The best fit model with flux 
constraints is still consistent with the ellipsoidal model. The 
fit statistics and the F-test results are shown again in Table [2] 
Note that the major axis of the mass is misaligned with the 
light when we use the flux ratios as a model constraint. 



5. SUMMARY 

We modeled the lensed quasar images and the Einstein ring 
formed from its host galaxy in the lens PG1 1 15+080 using a 
scale-free potential with arbitrary angular structure. The best 
fit ellipsoidal model is consistent with previous results, and 
provides a reasonable fit to the astrometric constraints of both 
the Einstein ring and the images. Non-ellipsoidal models are 
constructed by adding higher-order deviations from the ellip- 
soidal model, but none of the additional degrees of freedom 
improves the models. The best-fit non-ellipsoidal models are 
still consistent with the best fit ellipsoidal model. We also 
find that the center of mass of the lens is consistent with the 
measured center of light even when this is not imposed as a 
constraint. 

When we try to fit the fluxes as well, including the anoma- 
lous A1/A2 flux ratio, the best-fit non-ellipsoidal models still 
fail to match the observed flux ratios and are still consistent 
with the ellipsoidal model. We conclude that the suggestions 
that complex angular structure in the lens galaxy can explain 
the anomalous flux ratio in PG1 115+080 are incorrect - the 
galaxy is indistinguishable from an ellipsoid. This result does 
not address the problem of whether the anomaly is due to mi- 
crolensing or satellites and whether the source of substructure 
is in the lens or along the line of sight. The predicted A1/A2 
flux ratio for the astrometry only models is 0.97 ± 0.03 (see 
Tab[Q. There is some evide nce that the emission lin e flux 
ratios are closer to this value (Popovic & Chartas 2005), sug- 
gesting that the anomaly is due to microlensing rather than 
satellite despite the lack of time variations in the ratio. Since 
Einstein rings are relatively common in H-band HST obser- 
vations of lenses, we should be able to test the ellipsoidal hy- 
pothesis in many additional systems. 
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